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Abstract 

X: 

Using a program of perturbative resummation I compute the damping 
rates for fields at nonzero spatial momentum to leading order in weak cou- 
pling in hot QCD. Sum rules for spectral densities are used to simplify the 
calculations. For massless fields the damping rate has an apparent loga- 
rithmic divergence in the infrared limit, which is cut off by the screening of 
static magnetic fields ("magnetic mass"). This demonstrates how at high 
temperature even perturbative quantities are sensitive to nonperturbative 
phenomenon. 



I. Introduction 



For an asymptotically free theory such as QCD, at high temperature perturbation 
theory is a reasonable first approximation. Even if initially there are no bare masses, 
in an interacting plasma mass scales small relative to the temperature T are generated 
radiatively. With g the QCD coupling constant, elementary diagramatic techniques 
show that quasiparticles acquire thermal "masses" of order gT at one loop order. 
For example, both time dependent electric and magnetic fields are screened by a 
thermal gluon mass m g ~ gT, as are static electric fields. Static magnetic fields are 
not screened to this order, since the plasma is one of (colored) electric charge. It is 
expected that static magnetic fields are screened nonperturbatively by a "magnetic 
mass" 0. 

The thermal masses are related to the real part of the pole in the quasiparticle 
propagator. The imaginary part of the pole is proportional to the damping rate, 
7, and determines how rapidly a system near equilibrium approaches it. Unlike the 
screening lengths, which are easy to compute, even to lowest order in g the damping 
rates can only be computed consistently after the resummation of an infinite set of 
diagrams, termed hard thermal loops While infinite, the entire series of hard 

thermal loops can be succintly expressed in terms of simple effective actions . This 
resummation program has been applied to a variety of problems ||3|-[5], §)-|T^]- 

The damping rates are inevitably of order g 2 T, but they depend in an interesting 
fashion on how fast the quasiparticle is moving through the thermal medium [0. 
As a consequence of Landau damping, resummation produces an effective gluon that 
propagates below the light cone, with damping dominated by scattering off of spacelike 
effective gluons. If the incident field is initially at rest, both the recoil field and the 
effective gluon carry nonzero energy and momentum. While the calculations are 
involved, the effective gluon only probes energies and momenta of order gT, and so 
the damping rate is just some pure number times g 2 T ||,|],|]||. 

If the initial field is in motion, however, the effective gluon can be emitted at 
ninety degrees relative to the incident (and final) direction. When the field is moving 
sufficiently fast, the effective gluon carries almost zero energy and spatial momentum, 
and yet still contributes to the damping rate. In this way, transverse effective gluons 
probe the static magnetic sector, and so are sensitive to the presence of a magnetic 
mass. For example, in an SU(N) color gauge theory the damping rate for a gluon 
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moving with momenta of order T is shown in sec. IV to be 

_ g 2 NT ft m 2 g \ + imQ ^ \ . (1 ^ 

8tt \ \m 2 mag + 2m mag 7 t J ' '"/ 

for technical reasons this expression only holds in the limit where m mag ^> j t . 

Over large distances the behavior of static magnetic fields is controlled by a purely 
bosonic effective gauge theory in three dimensions, with a (dimensional) coupling 
constant = g 2 T. Assuming that bosonic QCD in three dimensions has a mass gap, 
the magnetic mass must be a number times g 2 T, m mag ~ g 2 T, up to a possible factor 



of yln(l/g) (T^j. Thus the damping rate in (|Llf) is of order ^ t ~ g 2 Tln(l/g 2 ). It 
was first shown in ref. that factors of ln{l/g 2 ) are generic to damping rates at 
nonzero velocity: here I consider more carefully how the logarithm is cut off, and try 
to evaluate the constant under the logarithm. Thus while (1.1) holds for m mag ^> 7$, 
unless m mag is some very large number times g 2 T, this limit probably does not apply: 
more likely, m mag < j t . I present this result to demonstrate that it is possible, at 
least in principle, to compute the constant under the logarithm. 

The formula in (1.1) does not apply to hot QED. In hot QED the behavior of 
static magnetic fields is determined by bosonic QED in three dimensions. This is 
a free theory, so = 0, and the damping rate for fast fermions, analogous to 

(1.1), is logarithmically divergent. (The damping rate for fast photons is not very 
interesting: it is finite and of order e 3 T.) I suspect that the damping rate for fast 
fermions in hot QED is finite and of order e 2 T/n(l/e), but this requires separate 
analysis beyond that presented here. Amusingly, as suggested originally in ref. [[|, 
the calculation in hot QED is more difficult than that for hot QCD, at least when 
m mag > 7* fT§ . 



Notice that the damping rate appears in the argument of the logarithm. Lebedev 



and Smilga [[IJJ first pointed out that it is necessary to include "f t self consistently, 
which is how it enters into the right hand side of (1.1). Their calculations indicated 
that 7t alone suffices to cut off the logarithmic divergence, even if m mag = 0. Recently, 
Baier, Nakkagawa, and Niegawa [16] argued that while the damping rate must be 
included self consistently, that if jt is determined from the position of the singularity 
in the propagator, then without the magnetic mass, by itself jt does not cut off 
the logarithmic divergence. My result in (1.1) accords with their arguments. This 
happens because the analytic structure of the propagator is rather complicated, with 
an unexpected branch cut appearing off the physical sheet, near the pole in the 
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quasiparticle propagator. The pole has an imaginary part 7^ the branch cut begins 
at a point that is seperated by an amount m mag from this pole. The restriction that 
m mag 3> It arises because it is much simpler to treat the case when the branch cut 
and the pole are far from each other than if they are close. 

Linde first argued that a magnetic mass renders the free energy sensitive to non- 
perturbative effects at four loop order Jl|. The above shows that for fast fields the 
damping rates are sensitive to such effects even at leading order. I confess that while 
I introduce the magnetic mass in a plausible fashion, it is at best a caricature of 
nonperturbative effects. Moreover, the appearance of the magnetic mass does not 
imply that calculations are fruitless: perhaps a marriage of perturbation theory (as 
in (1.1)) and lattice gauge theory (to determine m mag ) can be arranged. 

In sec. II I discuss I derive some necessary sum rules and therefy introduce the 
magnetic mass. In sec. Ill the damping rate for a slow, heavy fermion is computed. 
For kinematic reasons this damping rate does not probe very small momenta, of 
order g 2 T, and so is insensitive to the magnetic mass or to the details of analytic 
continuation. The damping rate for fast, massless quarks and gluons is computed in 
sec. IV. In sec. V the Ward identities are used to compute the leading logarithmic 
dependence in the damping rates of quarks and gluons traveling with momenta greater 
than order g 2 T. An appendix discusses how the the sum rules of sec. II can be used 
to compute the term of order g 3 in the free energy. 

For calculational ease all damping rates are computed in Coulomb gauge. I appeal 
to general proofs of gauge invariance in refs. || and [|I7[] to establish that this result 



is independent of the choice of gauge. Recently, Baier, Kunstatter, and Schiff [18 



observed that naive calculation in covariant gauges appear to violate these general 



proofs. The dilema was resolved by Rebhan [PH[, who demonstrated that an infrared 
regulator is required to treat the mass shell singularities which arise in covariant 
gauges; see, also, refs. ||, p0| , and fll|. I avoid these delicacies by sticking with 
Coulomb gauge, but note that explicit calculations in other gauges may well require 
the introduction of infrared regulators. 



II. Sum rules 

The conventions and notation of ref. [[| are followed. For an SU(N) gauge theory 
with Nf flavors of massless quarks in the fundamental representation, the effective 
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gluon mass induced by the thermal medium is 

In Coulomb gauge the only nonzero components of the gluon propagator *A tlv are 
*A 00 (K) = *A e (K) and *A'^(K) = (5 ij - titi) *A t (K). Here the gluon four mo- 
mentum is = (k°,k) and k = k/k; analytic continuation to real energies, with 
k° = iu, is implicit. When the momentum is soft, with uj and k of order m g , the 
effective plasmon and transverse gluon propagators are given by the tree term, plus 
the corresponding hard thermal loop fl2l] , |22 



ikP ~~ 

T 



*AJ\K) = k 2 - 3m 2 g Q l {"—) , (2.2) 



VW-^-l^tft^)-*^)-!). (2-3) 
The Q n are Legendre functions of the second kind. 

The propagator determines the spectral densities of the effective fields. For the 
transverse field, 

*A t (k°, k)= dr e ik T / duo k) (1 + n{u)) e^ T , (2.4) 



o 



where n(u) = l/(exp(u/T) — 1) is the Bose-Einstein statistical distribution function. 
The transverse spectral density is given by 

%(u, k) = Im *A t (-ioj + + , k)/n , (2.5) 

and is a sum of pole and cut terms, 

%(u,k) = Z t (k) Uu-El) + 5(u + E*j) + j3 t {u, k) -&{k 2 - uj 2 ) ; (2.6) 



is the step function, i?(ar) = 1 for x > 0, = for x < 0. The delta functions 
in (|2.6p represent the propagation of transverse gluons as quasiparticles with energy 
uj = El and residue Z t (k). The spectral density also includes the contribution of a cut 
below the light cone, |a>| < k, with spectral weight flt(uj, k). This cut is the result of 
Landau damping in a thermal distribution. The plasmon spectral density, pi(u,k), 
is defined similarly from *Ag(K), and determines the plasmon mass shell Ej~, residue 
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Ze(k), and the plasmon cut, Pe(u,k). Complete expressions for these quantities are 



given in refs. |2T| and |p3|| . For example, about zero momentum 

0) = (- IjL] * Pt (u,h) = (S(u - m g ) + 6(u + . (2.7) 



o 



That is, at rest a transverse gluon and the plasmon are degenerate in energy, E [ 
Eq = m g . The residue for the plasmon is a bit unusual — it is proportional to 
— 1/k 2 — but this is innocuous Jffi] . The mass shells for transverse gluons and for 
the plasmon split away from zero momentum, and are determined numerically as the 
solution of transcendental equations. 

For later purposes, note that about zero energy the contribution of the cut terms 
to the transverse and plasmon spectral densities are 

— Ityi 

7r k 2 — 3iim 2 ui / (4k) 



o^O - Im U n^, ,n^ > (2-8) 



and 

1 1 
BAu.k) ~ n — Im — — — — . (2.9) 

H K ' Jw ^° 7T k 2 + 3m 2 g + 3nm 2 g uji/(2k) v ; 

The imaginary terms in each propagator are both due to Landau damping, propor- 
tional to m 2 u/k as ui — ■> 0. The real terms are given by limits of the propagators at 
zero frequency. Thus for the plasmon term k 2 + 3m 2 enters in the denominator, with 
3m 2 the static electric mass squared. For the transverse term only k 2 appears, since 
at this order static magnetic fields are not screened. 



In sec. Ill only the expressions in ( p.8[) and (|2.9|) are required. In see's IV and V 



we need integrals of the spectral densities with respect to powers of uj. These integrals 
can be evaluated by means of sum rules p3| , ^0| , |2^| . 



The derivation of sum rules is an elementary exercise in complex analysis. While 
familiar, I go through several examples in order to emphasize the relevant physics. 
The essential point is to turn the integral over u into a contour integral in the plane 
of complex k°; to avoid confusion I relabel complex k° as z. For example, 

/+oo ^ If 
uj p t (u>,k) duj = f z A t (z,k) dz . (2-10) 

-oo 2-7T2 JC 

The contour C runs counter clockwise around the imaginary z axis. Since there are 
no intervening poles, the contour can then be deformed into a loop at infinity. For 
large z, the hard thermal loop 5U(z, k) falls off as m 2 /z 2 , and the effective propagator 
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behaves as the bare one, *A t (z,k) ~ 1/z 2 . Thus ( |2.10| ) is the same as for free field 
theory, 

to p t (u,k) du = 1. (2-11) 



At zero spatial momentum this sum rule is dominated by the pole terms in the spectral 
density, and is easily checked by using ( |2.7| ). 

The sum rule in ( |2.11| ) is familiar as a consequence of the equal time commutation 



rules. It is only valid to lowest order in g 2 , when the effective propagator includes 
just the hard thermal loop. For example, if the effective propagator included the full 
gluon self energy at one loop order, then the right hand side of fl2.ll ) is modified by 



terms for the (ultraviolet divergent) wave function renormalization constant of the 
gluon. This is distinct from the finite renormalization constant, Z t (k), above. 

Two other sum rules are needed for what follows. One is a relation for the plasmon 
spectral density: 

"+°° du * , , . 1 / 1 * A / , \ , 

— Pe{u,k) = -— f - A e (z,k)dz 

= -i + *A,((U) = -i + jr^. (2-12) 

The contour in the complex plane is now C — O, where O is a circle about the origin. 
The modification of contour is required because of the factor of 1/z in the integrand: 
this factor generates a pole in z whose contribution must be included. The result on 
the right hand side is a sum of two terms. The first results from deforming C into the 
circle at infinity. As for the transverse density, the contribution from the hard thermal 
loop vanishes at large z, and so the integral over C gives the same result as in free 
field theory, —1/k 2 . Secondly, there is the contribution from O; there the residue of 
the integrand at z = is just the value of the plasmon propagator at zero frequency, 
A^(0, k) = 1/ik 2 + 3m^). About zero momentum this sum rule is dominated by the 
pole terms, (|2.7| ). 

The last sum rule required is superficially similar to that for the plasmon density: 

— %(w,k) = - ^— I - *A t (z,k)dz 

oo 10 ITXI JC-O Z 

= *A t (0,k) ee 1 ■ (2.13) 

ft ' m mag 
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The contour at infinity does not contribute because *A t (z, k) falls off as 1/z 2 at large 
z. For the contour about the origin, O, the residue of the integrand is equal to the 
value of the transverse propagator at zero frequency, *A t (0, k). 

Using the effective propagator of ( |2.3| ), which includes just the hard thermal loop, 
( |2.13|) equals *A t (0, k) = 1/k 2 . Unlike the two previous sum rules, about zero spatial 
momentum (|2.13| ) is dominated not by the pole term, ( p.7|) , but by the cut term in 
the spectral density, ( |2.8| ) p3 [. 

In ( |2.13| ) I extend this relation, and introduce the magnetic mass, by defining the 
magnetic mass as the position of a presumed pole in the static transverse propagator, 
*A t (0, k) = l/(k 2 + m^g). This is merely a crude parametrization of the complicated 
physics which is responsible for the dynamical generation of a finite correlation length 
for static magnetic fields. As written, the magnetic mass represents the effects of a 
single glueball; surely there is a entire tower of glueball states, none of which need 
show up simply as a pole in the transverse propagator. A better approach would be 
to relate the quantities which enter into the damping rate to the vacuum expectation 
values of gauge invariant operators, which could then be computed by lattice gauge 
theory. At present this noble goal is beyond my means. 

In the calculations of the hard damping rate, kinematics typically restricts the 
integral over the gluon spectral densitites to lie below the light cone. For example, 
for the sum rule of (|2.13|) , 

ffc d u * , , N 1 2Z t (k) 

That is, the sum rule allows one to exchange an integral over the cut in the spectral 
density, p t (u,k), for a function of and Z t {k). 

To incorporate m magi the limiting form of the spectral density in ( |2.8| ) becomes 

k) I Im k 2 + m 2 mag -\nm 2 g u Ji /(Ak) ' (2 - 15) 

Using just (|2.15|) , the integral dup t (uj, k)/uj ~ l/(^ 2 + m ma g )- F° r & ~ ^mag 
this is the dominant term on the right hand side of ( [2.14P ; the pole term contributes 
~ 1/in 2 , which is smaller by order g 2 . This shows how the sum rule in (|2.13| ) is 
dominated by the cut term at small momenta. When k 3> m ma g the magnetic mass 
is negligible, and we recover the full sum rule. 



7 



III. Damping rate for a slow, heavy fermion 



For a heavy fermion of momentum P and mass M, the bare inverse propagator is 
Aj x (P) = — i P + M. To leading order the fermion self energy is 

E f {P) = - g 2 C f tr (7" *A f (P-K)Y *A^(K)) . (3.1) 

For a fermion in the fundamental representation the Casimir constant Cf = (N 2 — 
l)/(2iV); tr represents the integral over the loop four momentum K. In (|3.1| ) I have 
replaced the bare propagator by an effective propagator, *Af. This is defined as 
follows. For the bare propagator, the spectral density is 

p f (u,k) = (-0,7° + ik-j + M) ±.(6(u-E?) - 6(u> + E?)) , (3.2) 



where Eff = \/k 2 + M 2 is the fermion energy. To include the effects of damping I 
replace the sharp delta function in the spectral density by a Breit-Wigner form, with 
width 7j: 

(3.3) 

By the properties of the delta function this reduces to ( |3.2j ) as 7/ — > 0. The damping 
rate 7/ turns out to be of order g 2 T; this is g times the natural scale for the gluon 
spectral densities, which is set by the thermal gluon mass, m g ~ gT. This inclusion 
of higher order effects in a hard propagator extends the program of resummation 
outlined in ref. ||. I discuss later why it is valid to include these higher order 
effects, and not others, after computing 7/. Lebedev and Smilga [14j were the first 
to introduce the damping rate in this way. 

To simplify the computations I assume that the particle's motion is nonrelativistic, 
with a velocity u = p/M < 1. The case of relativistic motion is treated following the 
analysis of sec. IV. In ( pTT| ) the Saclay method is used to perform the sum over 
k°. The damping rate is proportional to the the imaginary part of the self energy on 
the mass shell. This is a sum of two terms, from the plasmon and transverse spectral 
densities: 

DiscE f {iE^,p) = 2«a t ( 7 ° - 1) - 10/(7° + 1), (3.4) 

where 

q 2 7rC f T f d 3 k f+°° du; f+°° , j f * , ,. „. u , 

(3.5) 
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Several approximations have been made to reach ( p.4|) and ( |3.5| ) . The term of leading 
order in g is given by replacing the Bose-Einstein statistical distribution function for 
the gluon by n(uj) ~ T/ui. Further, in most instances the fermion spectral parameter, 
oof, can be replaced by its average value, equal to the energy on the mass shell. At 
nonrelativistic velocities this energy is just the mass: oof ~ Ef_ k ~ M. The only 
instance where this is not allowed is in energy denominators: there, after taking the 
discontinuity of the self energy, one of the energy denominators produces the delta 
function for energy conservation in ( |3.5|) . (The other energy denominators don't 
contribute, since the spectral density for the soft gluon only has support for oo and k 
of order m g .) For that term the spectral parameter oof is redefined as oof — > 00f+E^_ k . 

Including the self energy, the renormalized fermion propagator is — i If + M — E/. 
The pole in the renormalized propagator is shifted from the bare mass shell to (Eff + 
ry/,p*), where 7/ is the damping rate, 

7/ = -2 a e + v 2 a t . (3.6) 

For simplicity, assume that the velocity, while small, is larger than g, 1 ^> v ^> g, so 
that Ef_ k ~ M + (p - k) 2 /(2M) ~ M - pkcosO/M. The delta function for energy 
conservation is used to fix the angle between p and k, cosO = p ■ k = (uj + Uf)/(vk), 

g 2 C f T r°° , 2 ,, r+°° duo r+°° , 7/ * . , , n . , . 
a* = ^ I kdkj^ - J^du,.— P i t{u ,k)Hvk-\ U + Uf \) . 

(3.7) 

The step function $ enters to ensure that \cos6\ < 1. I assume that the constraint 
is satisfied separately by \uo\ < vk and \uf\ < vk, so that the integrals over uof and 00 
decouple; this is justified following ( P-17| ). The integral over 00 f gives 

vk , If if vk\ 

duof 2 u 2 = 2 tan' 1 — . (3.8) 



■vk ' QOf+Jf \7/ 

The integral over the plasmon spectral density in ae is 



vk duo * . , . 3m 2 n v 

%(u,k) ~ - ——£—. (3.9) 



vk 00 (k 2 + 3rrig) 2 



Since vk in (|3.9| ) is small relative to k, the integral in ( |3.9|) is just 2vk times the limit 
of *pe{oo, fc)/w asu->0, ( |2.9| ). The remaining integral over k is finite, dominated by 
momenta k ~ m g ~ gT. Then I can take 7/ ~ in (|3.8j ), so that 

°« = - ht ■ < 310 > 
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This term is negative, and so from ( |3.6| ) a positive contribution to the damping rate. 
While various assumptions were made about the velocity to obtain ( |3.10|) , it is not 
difficult to go back to (|3.5| ) and show that one obtains identically the same result 
even for a field at rest, (5) of ref. 0. 

For the contribution of the transverse gluons, again since the integral runs only 
from — vk to v k, for small velocities only the limiting form in (|2.8| ) is required; the 
integral over it gives 

— %(u,k = tan- l [—-±\ . 3.11 

-vk to irk z \ 4k z J 

The magnetic mass, as enters in Q2.15Q , has been neglected; see the discussion fol- 
lowing (p.!7|) . Corrections to ( |3.11| ) are proportional to the velocity. After rescaling 
k — > jfk/v, ( |3.5|) becomes 

g 2 C f T r°° dk 1/t , i fc 2 \ . 

^ = hkl T tan ^ tan (*) • (312 > 

where c is a pure number, 

„ 3ii 2 m 2 v 3 

' 3 = -jj-- < 3 ' 13 ) 

In this instance it is necessary to keep 7/ 7^ 0, as otherwise ( pM2j ) develops a loga- 
rithmic divergence. Assume that the velocity lies in the range 

1 > v > g 2/3 . (3.14) 

Then the parameter c is large, since with m g ~ gT and 7/ ~ g 2 T, c 2 ~ m 2 v 3 / , ~f 2 ~ 
v 3 / g 2 ^> 1. It is then straightforward to compute the integral in (|3 . 1 2| ) . The integrand 
behaves like 1/k only for c> fc > 1 and is otherwise well behaved. Up to corrections 
of order ln(c)/c, 

a * = 9 -T^ l < c2 )- ( 3 - 15 ) 
Ibirv 

Altogether, Q, ( g3Dp , ( fTCJ ), and (£T5|) give 

At small velocities, inside the logarithm I can replace 7/ by its value at zero velocity 
to obtain 

g 2 C f T ( v, fl67t 3 (N + N f /2) v 3 \\ 
-" = ht[ l+ 2 ln { 3Cj ?JJ + -- (317) 
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The coefficient of the logarithm agrees with previous results 0, |16|| . 

I now justify separating the constraints on uj and ujf in ( |3.7| ). For the term 
involving the longitudinal spectral density there is no question; 7j can be sent to zero 
at the outset. For the transverse spectral density, the integral over uj is dominated 
by uj ~ vk; from the right hand side of (ft.ll|) , the relevant scale of momenta is 
k ~ y/vrrig, so the dominant frequences are uj ~ v 3 ^ 2 gT. For v ^> g 2 ^ 3 , then, oj 3> g 2 T, 
and so the scale for uj is much greater than that for ujf, which is ujf ~ 7/ ~ g 2 T. 
This separation in scales produces the logarithm in ( [3.16| ) and (13.171) , and justifies 
treating the constraints separately. This is not allowed if Q3.14| ) is not obeyed. For 
smaller velocities, v < g 2 ^ 3 , the scales in uj and ujf do mix; ultimately there will be 
no logarithm, with the term from the transverse spectral density vanishing smoothly 
as v —>■ 0. Larger velocities, v ^> 1, puts us in the relativistic regime, which is the 
subject of the next section. 

Also, since the dominant momenta for the transverse density are k ~ \fvm g 3> 
g x l 3 m g ~ m mag / g 2 ! 3 , it is permissible to neglect the magnetic mass, taking (|2T8| ) for 
the limiting form of the spectral density instead of ( |2.15| ). 

I have been somewhat careless on one other point. The correct damping rate 
is given by evaluating the imaginary part of the self energy at the position of the 
pole in the propagator. Including damping, this pole is off the physical sheet, at 
ujpoie = Ep 4 + ijf- Instead, I evaluated the imaginary part at uj = E^ 1 , and assumed 
that the continuation to uJ po ie is trivial. Baier, Nakkagawa, and Niegawa |16| have 
recently argued this continuation can produce a nonzero contribution to the damping 
rate. 

For velocities which satisfy ( |3.14|) , though, these subtleties can be overlooked. 
Suppose I were to evaluate the self energy not just for uj = E^ 1 , but for for uj = 
Ef + 6E, with 8E of order g 2 T. This alters energy conservation, so that in (|3.7f ) |u;| 
and \ujf \ must be < v k + 5E. For 5E ~ g 2 T, however, this change is negligible, since 
vk ~ v 3 / 2 gT g 2 T if v ^> g 2 ^ 3 . Thus the continuation from 5E = to 5E = ryj 
does not affect the result for the damping to leading order in g. 

Consequently, for a slow, heavy fermion, due to kinematic reasons there is no 
sensitivity to the magnetic mass (contrary to what I claimed in ref. [0) nor to details 
of analytic continuation (unlike ref. [|16j, at least for velocities as in (|3.14|) ). Both of 
these effects do enter for a field moving at relativistic velocities. 
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Why is it that corrections from the imaginary part of the self energy must be 
included, and not those from the real part? While there are certainly corrections to 
the mass shell of order g 2 T, at hard momenta these are independent of the spatial 
momentum. By energy conservation in ( |3.5| ), however, at this order all that enters 
into the damping rate is the difference in energies, E^ — E^_ k — so a constant shift in 
the mass shell cancels out. Similarly, consider the contribution of those higher loop 
diagrams which can be represented as vertex corrections to the diagram at lowest 
order. Assuming that the vertex correction is fixed by the appropriate Ward identity, 
if the self energies are slowly varying functions of momenta, the vertex corrections 
are even better behaved, and so will not generate the type of logarithmic divergences 
found above. (These arguments do not apply at soft momenta, where the self energies, 
and so the vertex corrections, depend nontrivially on momenta.) Nevertheless, it 
would be well worth checking these naive arguments by explicit calculation at two 
loop order. 

IV. Damping rates for fast fields 

Since it enters accompanied by an overall factor of the velocity, the logarithmic 
sensitivity found for the damping rate of a slow, heavy field is relatively innocuous. 
The damping rate of a fast particle is much more sensitive to small momenta. I first 
consider the case of a transverse gluon at "hard" momenta, p ~ T. 

For a transverse gluon the bare propagator is A i T 1 (-P) = (p ) 2 + p 2 - About the 
mass shell E* = p, the bare propagator behaves as A^ 1 ~ —Zf 1 (p)(uj — E£), where 
Z t (p) = l/{2p) is the residue for a hard field. Following Baier, Nakkagawa, and 
Niegawa ||16|| , the damping rate is determined from the discontinuity of the self energy 
at the position of the singularity in the propagator. To take this into account, I 
evaluate the discontinuity not just at, but near the mass shell, by introducing the 
function 

T t (5E) = Z t (p) DiscU t (i(E* + 6E),p) . (4.1) 

Here IT t is the transverse part of the gluon self energy, iP J '(P) = (5^ — p l fP)I\.t{P) + 
. . .. I first evaluate the function T t (5E) for real SE, and then analytically continue 
in SE to determine the damping rate from 7 t = r^ryt). Implicitly, this definition 
assumes that the propagator has a true pole at u = E^ + ij t , which I shall show is 
justified in hot QCD when m mag ^ 0. In hot QED, it is necessary to define the 
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damping rate by the position of the singularity in the propagator; the singularity will 
not just be a simple pole. 

By the resummation proceedure of ref. 0, to determine the damping rate at hard 
momentum requires the evaluation of the self energy which differs only slightly from 
the usual one loop diagram. Kinematics requires that only one line in the loop is 
soft, so bare vertices can be used. With bare vertices, the only contribution to the 
discontinuity is from the diagram with three gluon vertices. (The quark loop does 
not contribute to this order because there is no enhancement from Bose-Einstein 
statistics.) After projecting out the transverse piece of the gluon self energy, 

T t (6E) = g 2 NZi\p)tr ((1 - (k ■ p)f *A t (K) ■ *A e (K)) A t (P - K) . (4.2) 



For the soft field, with spectral parameter uj, the spectral densities are those of sec. 
II. For the hard field, with spectral parameter Uh, the bare spectral density is replaced 
by 

pM ' p ' k) = ^-^(h-vw + (<* + i»+J ' (4 - 3) 

After doing the sum over k° and retaining only terms of leading order in g, 

d 3 k r+°° du 



USE) = fNT /— 3 / - / 



duo 



(2tt) 3 J-oo id J-oo LUf + It 

((1 - (k ■ pf) %(u, k) - *p e (co, k)) 5(cu + co t -E t p -SE + El_ k ) . (4.4) 

The spectral parameter for the hard field has been shifted by u t — > uj t + E*_ k . With 
the mass shells for ultrarelativistic fields, the delta function for energy conservation 
fixes the angle between p and k as cosO = (u + u t — 5E) / k. 

Again I assume that the integrals over ui and ui t decouple. The integral over u t is 
peaked about u t ~ j t ~ g 2 T, so there the effects of SE ^ must be included. The 
integral over the longitudinal spectral density, as in ( |2.12| ), is perfectly well behaved, 
allowing me to set 7 t , SE, and m mag to zero; this is also valid for the integral over 
u)p t {uj,k), as in (|2.11|) . The integral over p t (u , k) / 'u , as in Q2.13D , is in principal 



sensitive to ujf, SE, and m ma9 when k ~ m mag . I now make the further assumption, 
however, that 

m rnag > 7 t . (4.5) 
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The limits of integration over u run properly over ±k — uo t + SE. When (|4.q ) holds, 
however, even for k ~ m ma9 I can neglect the effects of w t ~ 5E ~ 7^ and just let 
the integral over uj run from ±/c, as in the sum rule of Q2.14|) . 

The integral over u t is 

k+5E , It 1 fk + 5E\ l fk-8E\ . . 

du t „ - 9 = tan" 1 + tan' 1 . 4.6 

-k+SE Ut+lt V It V It J 



This is similar to (|3.8|) , except that here v = 1 and the shift in the mass shell, SE, 
enters. 

The integrals over the soft spectral functions are done using the sum rules of 
( [2. 11| ), ( |2.12|) , and (|2.13|) . As in (|2.14j ) these sum rules are used to trade an integral 
from ±k for terms which involve the right hand side of the sum rule and pole terms. 
In this way T t = Y s t mg + T r t e9 is written as a sum of a singular term, 

g 2 NT r°° , „ { 1 



V s ™ 9 (5 E) 



47T 



TT \ V It 

and a regular term, 




+ I ( tan- 1 ( *±*^ ] + tan' 1 I — — I I j^-^- I , (4.7) 



mag , 



* 4vr io VP + m^ ^ 

- 2Z lW ^) + + -gpj) = 1.0%81.... (4.8) 

So that each integral is finite at large momentum, a term proportional to kjik 2 +m 2 ) 
is subtracted from the integrand of the singular term, and added to that for the 
regular term. After doing so, each term is separately finite and well behaved for 
both small and large momentum. The singular term is sensitive to momenta of order 
k ~ g 2 T, so there the dependence on m mag , jt, and 5E must be retained. The regular 
term depends only upon momenta of order k ~ m g , and so up to corrections of order 
g, it is a pure number. This number, T r t eg ~ 1.09681...., was determined by numerical 
integration. 

The analytic form of the singular term was computed in the following manner. At 
zero magnetic mass it is easy to show that 

n 2 NT ( m 2 \ 
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Next, I compute the derivative of Y s t in9 with respect to the magnetic mass squared 
dVf n9 {8E) g 2 NT ( 7t \ /•+« dk 



dm 2 mag 8^ V ttJJ-oo (k* + ml ag )((k + 6Ey + ji) ' 



(4.10) 



This is a standard loop integral in one dimension, and can be done using the Feynman 
parametrization for the denominators. After doing the integrals over dk and the 
Feynman parameter, a relatively complicated form for dY^ 9 jdm 2 mag results. Doing 
the integral over m mag , and knowing Y s t m9 reduces to (f4.9|) for m mag = 0, gives a 
simple result, 

r? n9 (6E) = 9 -^Lln[-, ^ —) . (4.11) 



Having computed for real 5E, the analytic continuation of (|4.11| ) to complex SE 
is evident. There are branch points at 

5E = ±i( 7t + m mag ) , (4.12) 

which is off the physical sheet. The damping rate is evaluated at the pole in the 
propagator, at 5E = ij t . In the limit when m mag ^> 7 this pole is well seperated from 
the branch point, and the damping rate is just 

g 2 NT ( ( m 2 



lt = rr 9 (i lt ) + TT 9 = In - 9 + 1.09681... , 

x - 1 1 "' 2 ag + 2m maglt J J 

(4.13) 



u " \ \""mag 1 """mag It/ 

which is the result quoted in (|1.1|) . 

The same manipulations can be carried through for the damping rate of a (mass- 
less) quark field at hard momentum. I denote the damping rate of the standard mode, 
for which the chirality equals its helicity, as 7+. In the end the only change is in an 
overall factor for the Casimir of the representation: 



8tt \ \m 2 mag + 27+m ma9/ 



7+ = In — 9 £ + 1.09681... . (4.14) 



These damping rates are the only ones of significance at hard momenta. At low 
momenta the gauge field and the quark fields each have collective modes, the plasmon 
and plasmino, respectively. But the residues of these fields are exponentially small 
for hard momenta, so these fields can be neglected. 
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What of the damping rate of a fast fermion in hot QED, where the magnetic mass 
vanishes? In this instance the above approximations are inconsistent: the position of 
the branch point in the propagator, as in (4.10), coincides with the position of what is 
supposed to be a pole. Thus the spectral density for the hard fermion is not the pole 
of a Breit-Wigner form, but a branch point. Presumably the damping rate, defined 
as the imaginary part of the position of the branch point in the propagator, is gauge 
invariant and of order e 2 Tln(l/e). The question of gauge dependence is now more 
involved. For instance, at zero temperature it is known that the fermion propagator 
has a branch point singularity at the electron mass; while the position of the branch 
point is gauge invariant, the strength of the singularity is not |2l| . 

The analysis of hot QCD in what is probably the realistic case of m mag < j t is 
even more involved. Then the propagator has a pole at p + ijt, and a branch cut 
beginning at p+i(m mag + 'Yt)- The spectral density for the hard field must now include 
the effects of both the pole and the nearby branch cut. 



V. Damping rates for light fields 

The logarithm in the damping rate of a fast field arises from a very limited kine- 
matic region: in the one loop diagram, one line is very near the mass shell, while the 
other line carries almost zero momentum. In this section I analyze the same kinematic 
regime for fields moving at momenta comparable to the scale of the thermal mass. I 
just calculate the coefficient of the logarithm, which is relatively simple to compute. 

According to the effective expansion ||, for soft momenta the bare propagator 
is replaced by one which includes the hard thermal loop, with the effective plasmon 
and transverse propagators are those of ( |2.2| ) and ( |2.3| ). For instance, about the mass 
shell u = Ep, the effective transverse propagator behaves 

VC^p) ~ -Z^{p){u - E$, (5.1) 

where Z t (p) is the residue. The leading corrections to the damping rate are then 
determined by an effective gluon self energy, TP 1 ' , from which the transverse term 
*R t is extracted as usual. Using ( |5.1|) , the pole of the corrected transverse propagator, 
*Af 1 — *n t , then determines the damping rate for a soft transverse gluon, jtip) to be 

7t (p) ~ Z t {p) Disc*n t (iE*p) . (5.2) 
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This is similar to the function introduced in ( f4.1| ). Since I don't compute the constant 
under the logarithm, the subtleties of the previous section can be ignored, and it 
suffices to evaluate the discontinuity on the effective mass shell, without including 
the damping rate. The damping rates for the plasmon, and the quark modes, are 
defined similarly, and given at the end of this section. 

From (4.21)-(4.25) of ref. @, the effective gluon self energy is a sum of three 
terms, 

*tp"(p) = *n(£( p ) + + n ^( p ) » ( 5 - 3 ) 

where *n 3g is the graph with three gluon vertices, 

"W) = ^ftr soft *r"*(-P + K,P,-K) *A XX '(K) 

*T x ' ua '(-K, P, -P + K) *A a ' a (P - K) . (5.4) 
and *n 49 involves the four gluon vertex, 

*n%(P) = - y tr soft *r^(P, -P, K, -K) *A X °(K) . (5.5) 

These two terms are the same as in the bare expansion, except that bare propagators 
and vertices are everywhere replaced by effective quantities (the propagators are all 
in Coulomb gauge). Lastly, there is the ghost loop, IL^; because there are no hard 
thermal loops in ghost amplitudes, this loop equals that in the bare expansion. In 
Coulomb gauge the ghost loop has zero discontinuity, and so doesn't contribute to 
the damping rate. The subscripts on the trace in ( |5.4j ) and ( |5.5| ) indicate that the 
dominant term is given by the integral over soft momenta. 

The effective vertices which appear in *Yl are nontrivial functions of momenta, 
and so the discontinuity of such diagrams is far more complicated than in the bare 
expansion. For example, with a bare vertex the tadpole diagram has no discontinuity. 
In contrast, the discontinuity of *n.4 g is nonzero, because the vertex itself has a 
discontinuity from Landau damping. 

Assume that the leading logarithmic behavior of the damping rate at nonzero 
momentum arises from the same kinematic regime as for a fast field. Then while the 
diagram with a four gluon (effective) vertex, *Il4g, does contribute to the discontinu- 
ity, it can't generate a logarithm, since after cutting through one line and the vertex, 
there is no soft line left to integrate over. Thus the only diagram to contribute is that 
with three gluon vertices, *Il3 9 . If the gluon with momentum K is very soft, with 
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K ~ 0, then the other gluon, with momentum P — K, is very near its mass shell. (Of 
course I have to muliply by two, since the two gluons could be interchanged: K could 
be near P, so P — K is very soft.) Now much of what makes computation in the 
effective expansion so involved is the momentum dependence of the effective vertices. 
Under the assumption of very soft K, however, the Ward identities can be used to 
avoid having to compute any vertices whatsoever. Similar Ward identities have also 
been obtained by Weldon [B4[ . 



The effective three gluon vertex satisfies the Ward identity 

K x *T^\P,-P-K,K) = *A~'(P + K) - *A;J(P), (5.6) 
Hence as K — > 0, 

*r^(p,-p,o) = %1{P). (5.7) 

At small K (|5.7p can be used to compute the effective three gluon vertices which 
appear in *Tl3 9 , ( |5.4| ). Even more labor can be saved by recognizing that only part 
of 05.71) contributes. 



I introduce polarization vectors, e l a (p). These are defined as transverse to p, 

p'eOp) = 0, (5.8) 

and an or t honor mal set, 

4(p)4(p) = 5 ab . (5.9) 
They can be combined to form a projection operator in p as 

p*ffi = gt-tfpi = £ ei(p)ei(p). (5.10) 

o=l,2 

The advantage of using the polarization vectors is that by sandwiching the effective 
self energy between the e l 's, from ( |5.9| ) I automatically project onto the transverse 
part. Further, for the gluon of momentum P — K in the loop, in Coulomb gauge the 
propagator can be approximated as 

*A ij (P - K) = P i] (p- k) *A t {P -K)~ P ll {p) *A t (P - K) . (5.11) 

Then (|5.10|) is used to write P^(p) as a sum over polarization vectors e l a (p). 

In this way, all I need of the Ward identity in ( |5.7j ) are the terms which survive after 
sandwiching it between polarization operators. Using the behavior of the transverse 
propagator in (|5.1|) , 

e* a (p) *r^' fc (P, -P, 0) el(p) = 5 ab p k Z-\p) v t (p) . (5.12) 
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Here Vt(p) is the group velocity for a transverse gluon of momentum p on the effective 
mass shell, 

OE * 

v t (p) = (5.13) 

Of course at hard momenta Vt(p) — > 1. 

Computing in the kinematic regime when one transverse gluon is very soft, and the 
other gluon is transverse and very near its mass shell, it is then trivial to mimic the 
calculations of the previous sections to extract the leading logarithm in the damping 
rate. For the gluon with spectral parameter uj t , near its mass shell I replace the delta 
functions for the pole term by smeared Breit-Wigner forms with width j t . If the soft 
gluon has spectral parameter uo the contribution to the damping rate from a very soft 
transverse gluon is 

d 3 k r+°° duo 



, , 9/ , f d°k r +co duj f +co , 

lt {p) ~ g 2 NT v 2 (p) Tns / — / duo. 

Jk<t:m„ (Z7TI J-oo UJ J-oo 



It 

t , ,2 i „,2 



fc<m 9 (27r) 3 J-oo UJ J-oo UJ 2 + ^ 

[l-(k- Pf ) %(uj, k)) 5(uj + uh- El + El_ k ) . (5.14) 

This is almost the same integral as for the damping rate at hard momentum ( |4.4j ), 
except that each vertex contributes a factor of the group velocity, and the mass shell 
is now that for an effective field. Expanding the delta function for energy conservation 
in small k fixes the angle between p and k to be cosO = (uj + uj t ) / (vt(p)k): notice the 
extra factor of one over the group velocity, which is like the nonrelativistic case in 
sec. III. Using the delta function to integrate over 6 gives one factor of l/\v t {p)\, so 
that in all 

g 2 NT ( 1 \ 

lt,t{P) ~ Vt '^ P ' ln \~g 2 ) + '" ' 

I approximate the argument of the logarithm as ^n(m 2 /m 2 nag ) ~ ln(l/g 2 ). As indi- 
cated in ( |5.15| ), the damping rate of the plasmon can be computed similarly: it is 



given by replacing v t {p) with the plasmon group velocity, vg(p) = dE^/dp. 

For the quark field there are two modes: at positive energy, the standard mode 
has chirality equal to helicity, and while for the plasmino, its chirality is equal to 
minus its helicity ||26|| . Denoting the mass shells by E^ 1 , in analogy to the result at 



hard momenta, ( |4.14| ), the damping rate of the quark field is 



~.±M>) - ^-^Mp)|/nl^) h ... . (0.10) 
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where v±(p) = dE^ /dp are the group velocities for the quark modes. In this instance 
I explicitly write the absolute value of the group velocity. This was not necessary 
before, since the mass shells for the transverse, plasmon, and standard quark modes 
are each monotonically increasing in p, so the group velocities are always positive. The 
plasmino mass shell, however, decreases from zero momentum, reaches a minimum at 
p = p c , and then increases, so is negative for p c > p > 0. Similarly, 7_(p c ) in 

( p,16| ) vanishes at p = p c since v_{p c ) = 0. 

The terms in ( 5.15 ) and ( |5.16|) apply for momenta p ^> g 2 T. For momenta p ~ T, 
all group velocities v — > 1, and the results of the previous section are recovered. The 
restriction that the momenta be greater than g 2 T is necessary because I assumed 
that I could approximate E p — E p _j, ~ v(p)cos8, which is incorrect for p <C m g . For 
instance, if one computes as above, but now at exactly zero momentum, p = 0, one 
finds no term as ln(l/g 2 ) in 7(p). The crossover scale at which a ln(l/g 2 ) appears in 
the damping rate is set by the width of the spectral density for a field near its mass 
shell, which is of order ~ g 2 T. This restriction is of no concern for the gluon fields, 
since v tt e(p) ~ p/m g as p — > 0. For the quark fields, however, |f±(p)| — > 1/3 as p — > 0, 
and this caveat is important: there is no ln(l/g 2 ) in 7±(0), as indicated by the naive 
extrapolation of ( |5.16p . 

This is consistent with explicit calculations at zero momentum. The damping 
rates of both gluons 0] and quarks at zero momentum have been computed: 
both results are a pure number times g 2 T, with no terms as ln(l/g 2 ). Similarly, the 
damping rate of the plasmino at the minimum in its dispersion relation, p = p c , is 
surely a pure number times g 2 T. 

It would take some effort to compute the constants under the logarithm in the 
damping rates of ( |5.15| ) and ( |5.16| ). Not only does the other diagram, *n 4s , enter, 
but for both diagrams the full form of the vertices are required. Also, as in sec. IV, 
the effects of analytic continuation to the pole in the propagator must be included. 



Appendix: The plasmon pressure and sum rules 



van Weert and collaborators [^7] raised the following question: the order g 3 term 
in free energy arises entirely from the plasmon through the static electric mass. The 
transverse modes do not contribute because static transverse fields are not screened 
perturbatively. On the other hand, for time dependent fields both the plasmon and 
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transverse fields are screened by thermal masses m g . So why don't transverse fields 
fields contribute to the free energy at order g 3 l In accord with an analysis by Toimela 
28j , in this appendix I show how the sum rules can be used to demonstrate the 



cancellation of the transverse pressure at order g 3 . 

The terms of order g 3 in the free energy arise from the summation of "ring" 
diagrams [29|, as the free energy of effective fields. For the transverse gluon, per color 
degree of freedom the two transverse modes give a free energy 



% = tr (inCA^iK)) - ln(K 2 )) ; (a.l) 



as before, tr denotes the integral over the four momentum K . The effective propagator 
*A^ 1 (K) = K 2 — SHt(K), where 5Ut(K) is the hard thermal loop in the transverse 
gluon self energy. From ( |2.1| ) and ( p.3|) , 5U t (K) is g 2 T 2 times a function of momentum. 
Then it is easy to compute the derivative of *F t with respect to g 2 : 

The sum over k° is done by using the spectral representation of the effective propa- 
gator, 

/ \ r d 3 h /"+0O 

tr(K 2 *A t (K)) = J-^- J ^ duj (1 + n(uj)) (-uj 2 + k 2 )%(uj,k) . (a.3) 

The dominant term is given by approximating the Bose-Einstein statistical distribu- 
tion function n(uS) ~ T/u, 

tr %M) ~ T J H /_ + ; g k) . (,4) 

Now use the sum rules of ( |2.11| ) and (|2.13| ) to do the uj integrals, 

Hence if the magnetic mass is ignored, the contribution of transverse gluons to the 
pressure cancels identically. With m mag ~ g 2 T, up to a possible Jln(l/g 2 ) the 
first nonvanishing contribution to the transverse free energy is of order *F t ~ g 6 T A . 
This is where the effects of the magnetic mass are expected to arise in the free energy 
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Why did this cancellation occur? The arguments of van Weert et al |^7] considered 
only the contributions of pole terms to the spectral density. This is the dominant 
contribution to the sum rule of ( |2.11|) . But the effective spectral density of (|2.6| ) also 
includes cuts from Landau damping, and it is these which dominate the sum rule of 



( gig ). For 



■may 



in *F t , the contribution of the cuts identically cancels that of 



the pole terms. 

The analogous contribution of the plasmon to the free energy is, per color degree 
of freedom, 



1 



tr 



(Zn(*A7 1 (K)) - ln(k 2 



(a.6) 



where *A e 1 (K) = k 2 — £ILj (K), with SXl^K) the hard thermal loop in the plasmon 
part of the gluon self energy, Q2.2| ). The derivative of this term with respect to g 2 is 



d% 
dg 2 



2<r 



tr(l - k 1 *A e (K) 



Then 



Using the sum rule of (|2.12|) 



T 
2 



(2k) 



d 3 k 



-°° duo 

X> UO 



%{uj, k) . 



ln{k 2 + 3 m 2 



(a.7) 



(a.8) 



(a.9) 



Remembering that the static electric mass squared m el = 3m 2 , ( |a.9|) is exactly equal 
to the free energy as computed in the imaginary time formalism, including only the 
effects of the term with k° = 0. The result for ^ is a term of order g 2 (this is part of 
the free energy at two loop order) plus the term of interest, = — Tm? e J (12k) ~ g 3 T 4 . 
Thus sum rules demonstrate the equivalence with the results of the imaginary time 
formalism. 

The work of sec. II, and part of sec. IV, was done in collaboration with E. Braaten, 
to whom I am indebted. I also thank R. Baier and A. Weldon for discussions. This 
work was supported in part by the U.S. Department of Energy under contract DE- 
AC02-76CH00016. 
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